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We explore the viability of using solely spacelike information about a Schwinger function to extract properties 
of bound states. In a concrete example it is not possible to determine properties of states with masses > 1.2 GeV. 
Modern Dyson-Schwinger equation methods supply a well-constrained tool that provide access to hadron masses 
and cr-terms. We report values of the latter for a range of hadrons. Of interest is analysis relating to a u,d scalar 
meson, which is compatible with a picture of the lightest 0 ++ as a bound state of a dressed-quark and -antiquark 
supplemented by a material pion cloud. A constituent-quark <r-term is defined, which affords a means for assessing 
the flavour-dependence of dynamical chiral symmetry breaking. 


1. Introduction 

The pseudoscalar spectrum [ I G (J PC ) = 
1“(0 _+ )] contains three states below 2 GeV: 
7 t(140); 7t( 1300); and 7 t( 1800). The lightest has 
been much studied. However, a comprehensive 
understanding of QCD requires an approach that 
admits the simultaneous study of the heavier 
pseudoscalars and, indeed, other systems. An un¬ 
derstanding of the hadron spectrum and its real¬ 
isation within QCD is necessary in order to un¬ 
ravel the nature of the long-range force between 
light quarks. 

The lightest pseudoscalar meson is both a 
bound state of u- and d- quarks and the Goldstone 
mode in QCD associated with the dynamical 
breaking of chiral symmetry. This is readily and 
clearly understood using the Dyson-Schwinger 
equations (DSEs) [1,2]. The DSEs have proven 
a particularly useful device for studying the spec¬ 
trum and properties of light-quark systems. Mod¬ 
ern applications are reviewed in Refs. [3,4,5]. 

Herein we present results obtained from studies 
of the inhomogeneous Bethe-Salpeter equations 
(BSEs) for the scalar and pseudoscalar vertices 
in QCD. This is a practical means of mapping 
out the domain of applicability of the leading or¬ 
der term in a systematic and symmetry preserv¬ 
ing DSE truncation scheme [6,7,8]. Furthermore, 


we explore the capacity of such studies to com¬ 
plement contemporary numerical simulations of 
lattice-regularised QCD. 


2. Bound states from spacelike data 


The BSE provides a Poincare covariant tool 
with which to calculate the properties of bound 
states in quantum field theory. The inho¬ 
mogeneous equation for a pseudoscalar quark- 
antiquark vertex is 1 


r i(k-,P) 


n u = Z ^Y 


r A „ 


xi(q;P) K*™(q, k; P). 


( 1 ) 


where k is the relative and P the total momen¬ 
tum of the constituents; r,..., u represent colour, 
Dirac and flavour matrix indices; 


xi(q-,P) = S(q+)T :> 5 (q-P)S(q-) , (2) 

q± = q ± P/2; and / A represents a Poincare in¬ 
variant regularisation of the integral, with A the 
regularisation mass-scale [1,2]. In Eq. (1), K is 

1 For simplicity we work with two degenerate flavours of 
quarks. Hence, Pauli matrices are sufficient to represent 
the flavour structure. We employ a Euclidean metric with 
the conventions described, e.g., in Sec. 2.1 of Ref. [3]. 
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the fully amputated and renormalised dressed- 
quark-antiquark scattering kernel and S is the 
renormalised dressed-quark propagator. ( SSK is 
a renormalisation group invariant). The dressed- 
quark propagator has the form 

W _1 = z( p 2,( 2 ) P + M ^ ’ (3) 

and is obtained as the solution of QCD’s gap 
equation: 

Sip) -1 = Z 2 (*7 'P+ TObm) + S(p), (4) 

S(p) = Z\ J g 2 D liv (p-q)^-^ li S{q)Tl(q-,p) , (5) 
augmented by the renormalisation condition 

S(p)~ 1 \ p2=C 2=i'yP + m(C), ( 6 ) 

where m(() is the running current-quark mass at 
the renormalisation point (. These equations in¬ 
volve the quark-gluon-vertex, quark wave func¬ 
tion and Lagrangian mass renormalisation con¬ 
stants, Zi i 2 ) 4 (C, A )j each of which depends on the 
gauge parameter, the renormalisation point and 
the regularisation mass-scale. 

The solution of Eq. (1) has the form 

ir J 5 (k ; P) = y 7 5 [*£ S (fc; P)+ 7 • PF 5 (k ; P) 

+ 7 -kk- PG 5 (k-P) + a^k^P,, H 5 (k;P)] . 

( 7 ) 

This is the minimal and complete form required 
by Poincare covariance. The homogeneous pseu¬ 
doscalar BSE is obtained from Eq. (1) merely by 
omitting the driving term ; viz., Z. 475 2 ^-. The 
equation thus obtained defines an eigenvalue 
problem, with the bound state’s mass-squared be¬ 
ing the eigenvalue and its Bethe-Salpeter ampli¬ 
tude, the eigenvector. As such, the equation only 
has solutions at isolated timelike values of P 2 . 
On the other hand, the solution of the inhomo¬ 
geneous equation, Eq. (1), exists for all values of 
P 2 , timelike and spacelike, with each bound state 
exhibited as a pole. This is illustrated by Fig. 1, 
wherein the solution is seen to evolve smoothly 
with P 2 and the pole associated with the pseu¬ 
doscalar ground state is abundantly clear. 



Figure 1. Pseudoscalar amplitude Ps(0;P 2 ) 
obtained by solving the inhomogeneous BSE, 
Eq. (1), using the renormalisation-group- 
improved rainbow-ladder truncation described 
in [9]. The vertical dotted line indicates the 
position of the ground state tt mass pole. The 
inset shows l/Ps(0;P 2 ). 


A numerical determination of the precise loca¬ 
tion of the first pole in E^(k 2 = 0; P 2 ) will gener¬ 
ally be difficult. The task becomes harder if one 
seeks to obtain the positions of excited states in 
addition. It is for these reasons that the homoge¬ 
neous equation is usually used. However, if one is 
employing a framework that can only provide in¬ 
formation which is equivalent to the form of this 
vertex at spacelike momenta, then a scheme must 
be devised that will yield the pole positions. 

One obvious alternative is to focus on 
Pb(P 2 ) = 1 /E{k 2 = 0;P 2 ) and locate its zeros, 
and it is plain from Fig. 1 that this approach can 
at least be successful for the ground state. It 
is important to determine whether this (or an¬ 
other) approach can also be used in practice to 
determine some properties of excited states when 
information is available only for spacelike P 2 . 

While the DSEs can be used to generate such 
information, and the analysis of that information 
is what brought us to this point, herein for this 
purpose we consider a simple model for an in¬ 
homogeneous vertex whose analytic structure is 
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Table 1 

Parameters characterising our vertex Ansatz , 
Eq. (8). They were chosen without prejudice, sub¬ 
ject to the constraint in quantum field theory that 
residues of poles in a three-point Schwinger func¬ 
tion must alternate in sign [10], and ordered such 
that nii < rrii+i. We use b = 0.78. This is the 
calculated value of Z 4 (£ = 19 GeV, A = 200 GeV) 
used to obtain the curves in Fig. 1. 


i 

Mass 

Residue 

1 

0.14 

4.23 

2 

1.06 

-5.6 

3 

1.72 

3.82 

4 

2.05 

-3.45 

5 

2.2 

2.8 


known precisely; namely, 

M 

V(P 2 ) = 6 + ^-^- 
' P z + m 

i— 1 


( 8 ) 


where: for each i, rm is the bound state’s mass 
and a,i is the residue of the bound state pole in 
the vertex, which is related to the state’s de¬ 
cay constant; and b is a constant that represents 
the perturbative background that is necessarily 
present in the ultraviolet. The particular pa¬ 
rameter values we employ are listed in Tablet. 
This Ansatz provides a data sample that captures 
the essential qualitative features of true DSE so¬ 
lutions for colour-singlet three-point Schwinger 
functions, such as that depicted in Fig. 1. 

To proceed, we employ a diagonal Pade approx- 
imant of order N to analyse the data sample gen¬ 
erated by Eq. (8); viz., we use 


In(P 2 ) 


Co + CiP 2 + . . . + CnP 2N 
1 + C N+1 P 2 + . - . + C2nP 2N 


( 9 ) 


as a means by which to fit 1 /V(P 2 ). The known 
ultraviolet behaviour of the vertex requires that 
we use a diagonal approximant. NB. A real-world 
data sample will exhibit logarithmic evolution be¬ 
yond our renormalisation point. No simple Pade 
approximant can recover that. However, this is 
not a problem in practical applications because 
the approximant is never applied on that domain. 

In a confining theory it is likely that a colour- 
singlet three-point function exhibits a countable 


infinity of bound state poles. Therefore no finite 
order approximant can be expected to recover all 
the information contained in that function. Our 
vertex model exhibits M bound state poles. We 
expect that an approximant of order N < M will 
at most provide reliable information about the 
first N — 1 bound states, with the position and 
residue associated with the N th pole providing 
impure information that represents a mixture of 
the remaining M — (N — 1) signals and the con¬ 
tinuum. We anticipate that this is the pattern of 
behaviour that will be observed with any rank-IV 
approximation to a true Schwinger function. To 
explore this aspect of the procedure we studied 
the iV-dependence of the Pade fit. 

The domain of spacelike momenta for which in¬ 
formation is available may also affect the reliabil¬ 
ity of bound state parameters extracted via the 
fitting procedure. We analysed this possibility by 
fitting Eq. (9) to our Ansatz data on a domain 
(0,P 2 ax ], and studying the P 2 lax -dependence of 
the fit parameters. 

In all cases we found that a Pade approximant 
fitted to 1 /V(P 2 ) can accurately recover the pole 
residues and locations associated with the ground 
and first excited states. However, there is no 
cause for complacency. 

In Fig. 2 we exhibit the P^ ax -dependence of the 
mass-parameters determined via a N = 3 Pade 
approximant. Plateaux appear for three isolated 
zeros, which is the maximum number possible, 
and the masses these zeros define agree very well 
with the three lightest values in Table 1. This 
appears to suggest that the procedure has per¬ 
formed better than we anticipated. However, that 
inference is seen to be false in Fig. 3, which depicts 
the P^ ax -dependence of the pole residues. While 
the results for aq ,2 are correct, the result inferred 
from the plateau for 03 is incorrect. It is impor¬ 
tant to appreciate that if we had not known the 
value of 03 a priori, then we would very likely 
have been misled by the appearance of a plateau 
and produced an erroneous prediction from the fit 
to numerical data. Plainly, a N = 3 approximant 
can at most only provide reliable information for 
the first M — N = 2 bound states. 

We explored this further and applied aAf = 4 
approximant to the same model. In this case we 
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Figure 2. Pole positions (mass values) obtained 
through a fit of Eq. (9) with N = 3 to data for 
l/V(P 2 ) generated from Eq. (8) with the param¬ 
eters listed in Table 1. The coordinate ax is 
described in the text. Horizontal dotted lines in¬ 
dicate the three lightest masses in Table 1. The 
ground state mass (solid line) obtained from the 
Pade approximant lies exactly on top of the dot¬ 
ted line representing the true value. 

could still only extract reliable information for 
the first two poles. We subsequently biased the 
fit procedure by hard-wiring in Eq. (9) the residue 
and position of the lightest pole masses. This did 
not help. The N = 4 approximant could still not 
provide results that improved upon what we had 
already learnt with the N = 3 approximant. 

The results described herein, and our contin¬ 
uing analysis and exploration of other methods, 
including those popular in contemporary simula¬ 
tions of lattice-regularised QCD [11], suggest that 
solely from spacelike data it is only ever possi¬ 
ble to extract reliable information about bound 
states with masses which do not much exceed 
1 GeV. While these results are preliminary, they 
nevertheless provide sound reasons for caution. 

Following this background work on the viabil¬ 
ity in principle of using spacelike data alone to 
extract bound state information, we applied the 
method to the true DSE-calc.ulated pseudoscalar 
vertex that is in depicted Fig. 1. The decay con¬ 
stants and masses for the ground state pseu¬ 
doscalar and first radial excitation were obtained 
from the homogeneous BSE in Ref. [10]. The com- 


Figure 3. Pole residues obtained through a fit 
of Eq. (9) with N = 3 to data for 1 /V(P 2 ) gen¬ 
erated from Eq. (8) with the parameters in Table 
1. Horizontal dotted lines indicate the residues 
associated with the three lightest masses in Table 
1. The residue associated with the ground state 
(solid line) lies exactly atop the dotted line repre¬ 
senting the true value. The residue for the second 
pole exhibits a plateau at the correct (negative) 
value. However, the plateau exhibited by the re¬ 
sult for the residue of the third pole is wrong. 

parison between these masses and those inferred 
from the Pade approximant is presented in Fig. 4. 
As suggested by our background work, with per¬ 
fect (effectively noiseless) spacelike data at hand, 
reliable information on the first two states in this 
channel can be obtained. 2 We do not present a 
plot of the residues but they are accurate. In par¬ 
ticular, the ground state residue is positive and 
that of the first excited state is negative, and this 
is obtained without any bias in the fit. 


3. Sigma terms 


The cr-term for a state O is given by 


CTO 


= m( C) 


dmo 

dHO ’ 


( 10 ) 


where mo is the mass of the state, and it is a 
keen probe of the impact of explicit chiral sym¬ 
metry breaking on a hadron’s mass. The nu¬ 
cleon’s cr-term, <tn, has been estimated using: 

2 We are currently exploring the impact of Gaussian noise 
in the data. 
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Figure 4. Pole positions (mass values) obtained 
through a fit of Eq. (9) with N = 3 to the DSE re¬ 
sult for 1 /E 5 (k 2 = 0, P 2 ) depicted in Fig. 1. Hor¬ 
izontal dotted lines indicate the masses obtained 
for the ground and first excited state via a direct 
solution of the homogeneous BSE [10]. 


chiral effective theory, e.g. Refs. [12,13]; lattice- 
QCD, e.g. Ref. [14,15]; and QCD-based models, 
e.g. Ref. [16]. Our recent interest in hadron a- 
terms is motivated by their utility in using obser¬ 
vational data to place constraints on the variation 
of nature’s fundamental parameters [17]. In Ta¬ 
ble 2 we reproduce results calculated in Ref. [18] 
along with new results described in the text. We 
note that 

5mo = ero 5m(() 
mo mo m(C) 

so that the dimensionless quantity tabulated mea¬ 
sures the linear relative response of the mass of in¬ 
terest to a fractional change in the current-quark 
mass, which is a renormalisation group invariant. 

Radial Excitation of the Pion. We have a frame¬ 
work which enables the calculation of cr 7ri us¬ 
ing Eq. (10), where iri denotes the first pseu¬ 
doscalar radial excitation. To be specific, we can 
straightforwardly obtain the current-quark-mass- 
dependence of the 7iq in the renormalisation- 
group-improved rainbow-ladder (RL) DSE trun¬ 
cation described in Ref. [9]. In this truncation 


Table 2 

Calculated cr-terms. Those for n, N and A were 
reported in Ref. [18]. The values for p and u> 
listed herein were obtained via a direct analy¬ 
sis of the m(C)-dependence of the vector meson 
mass obtained in a solution of the rainbow-ladder 
truncation of the quark DSE and homogeneous 
meson BSE. They improve the values reported 
in Ref. [18], which were obtained using a simple 
fit to TO p (m(<C)) provided in Ref. [19]. All re¬ 
sults are renormalisation-point-independent and 
were obtained with: m U: d(() = 3.7 MeV, m s (Q = 
82 MeV, m c (C) = 0.97 GeV and m b (<) = 4.1 GeV, 
where C = 19 GeV. Perturbative evolution can be 
used to determine the associated renormalisation- 
point-independent current-quark-masses. 


H 

7r 

7Tl 

a 


<JH 

mu 

0.498 

0.017 

0.013 


H 

P 

to 

N 

A 

(Th 

m H 

0.021 

0.034 

0.064 

0.041 

q 


s 

c 

b 

<Jg 

Mg 

0.023 

0.230 

0.637 

0.851 


= 1.06 GeV [10] and we find 3 

off = 0.018 GeV. (12) 

This value is considerably smaller than that for 
the ground state pion: a n = 0.069 GeV. However, 
that merely serves again to emphasise the partic¬ 
ular character of the mass of QCD’s Goldstone 
mode and its amplification by dynamical chiral 
symmetry breaking. 

Scalar Meson . We do not pretend to have a 
complete understanding of the lowest mass 0 ++ 
state in the hadron spectrum. In this channel 
the rainbow-ladder truncation may be unreliable 
because the cancellations between higher-order 
terms in the systematic DSE truncation, which 
are so effective and important in pseudoscalar and 
vector channels, do not occur [20]. This is en¬ 
tangled with the phenomenological difficulties en- 

3 We are currently unable to provide a reliable estimate of 
meson-loop corrections to the mass and cr-term of the 7ri. 
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countered in understanding the scalar states be¬ 
low 1.4 GeV (see, e.g., Refs. [21]). 

Nevertheless, the rainbow-ladder DSE trunca¬ 
tion provides a light-quark scalar-meson solution. 
Combining the results from a number of sources, 
one finds [3] ro* L = 0.64 ± 0.06 GeV. With the 
kernel described in Ref. [9], which we have used 
throughout, one obtains 


to* l = 0.675 GeV, 

(13) 

= (0.184 GeV) 2 

(14) 

=> cr^ L = 0.025 GeV. 

(15) 


The scalar meson described by the rainbow- 
ladder truncation has a large coupling to two pi- 
ons [22]. It is therefore important to consider the 
effect of a 7T7 t loop correction to this state’s mass 
and cr-term. Such meson loop corrections can be 
estimated [14,18,23] and have a modest quanti¬ 
tative impact (< 15%) on <r p , ajy and cta and a 
larger effect (< 30%) on cr w . Their effect in this 
case can be analysed in the same way. 

We consider a single 7T7r-loop self-energy 


n 7 (ml) 


-3 9, 


2 

f77T7r 


16-7T 2 


dk 


k 2 u A „(k) 2 


uj(k) (w(fc) 2 — to 2 /4) 


,(16) 


where ui(k) = y/ m 2 + k 2 . The self energy both 
corrects m^ L and provides for the er —> 7T7 t width. 
In Eq. (16) 


3?nj 7r ((m^' L ) 2 ) = —(0.395 GeV) 2 . Hence, the 7T7T- 
loop acts to reduce the mass of the rainbow-ladder 
cr-meson and, from the shifted pole position, 

TO RL+7T7T _ q. 624 GeV. (19) 


This is an 8% reduction cf. Eq. (13). With all else 
kept fixed in Eq. (16), a variation of by ±20% 
alters the result in Eq. (19) by ±6%. NB. From 
Ref. [24], one might consider 0.60 ± 0.25 GeV as 
typifying the mass of the lightest scalar meson. 

We define the 7T7r-loop correction to the scalar 
cr-term via 


d 

dm 2 




(?7l^ ) ex pt. 


(0.161 GeV) 2 , (20) 


and combine this with Eq. (14) accordingly: 


m<c) aS(cT = <°- 090GeV > 2 

= ■ 2m^ 1_Ioop ^ a^ 1 " 100 ?) 

=> (t^ 1—loop ) = 0.0073 GeV, 


( 21 ) 

( 22 ) 


since to*° op ^ = 0.547 GeV. Equation (22) is a 
71% reduction cf. Eq. (15). The effect is large 
because the scalar-meson is very broad. The 
Y\y self-energy contributes much less to proper¬ 
ties of the rainbow-ladder p-meson because the 
width/mass ratio is significantly smaller [18,28]. 

It is noteworthy that a self-consistent solution 
of s — (to^ l ) 2 — nj 7r (3ffs) = 0 gives a pole position 


= 0.578 -i 0.311 GeV. 


(23) 


= 5.5i = 0.300GeV. (17) 

m“ L 

This is a typical value for the imaginary part of 
the T-matrix pole [24] and corresponds to 

T an ^ = 0.55 GeV, (18) 

which equates numerically to 0.92 (2 3y / s T ). 

The function u\ a {k 2 ) = 1/(1 + k 2 / A 2 ) 2 in 
Eq. (16) is a form factor, introduced to represent 
the nonpointlike nature of the it and cr and hence 
the cr7T7r vertex. The analysis of Refs. [25,26] 
indicates that the intrinsic size of the a de¬ 
scribed herein is 84% of that of the p. We there¬ 
fore choose a value of the regularisation mass- 
scale A ct = A pir7r /0.84. With A p7r7r = 1.23 GeV 
at m,r = 0.14 GeV, as determined in Ref. [27], 


The third iteration is the last to introduce a 
change >1%, and all iterations beyond the first 
reduce the real part in Eq. (23) by a total of < 7%. 
The effects are greater than that of the nucleon’s 
ttN self energy [18,29]. These notes provide a 
gauge for the accuracy of the one loop analysis. 

Constituent quarks. One measure of the impor¬ 
tance of dynamical chiral symmetry breaking to 
the dressed-quark mass function is the magni¬ 
tude, relative to the current-quark mass, of the 
Euclidean constituent-quark mass; viz., 

(. M E ) 2 := s 3 s = M(s) 2 . (24) 


For the current-quark-masses in Table 2 


Q 

u,d s c b 

M|(GeV) 

0.42 0.56 1.57 4.68 


( 25 ) 












7 



Figure 5. Ratio cjq/Mq [Eqs. (24) & (26)]. It is 
a measure of the current-quark-mass-dependence 
of dynamical chiral symmetry breaking. The ver¬ 
tical dotted lines correspond to the u = d, s, c 
and b current-quark masses listed in Table 2. 

A constituent-quark cr-term can subsequently 
be defined [18] 

dM E 

■ <26) 

It is a renormalisation-group-invariant that can 
be determined from solutions of the gap equation, 
Eq. (4). Our results are presented in Table 2. 

In Fig. 5 we depict <tq/Mq. It is a measure of 
the effect on the dressed-quark mass-function of 
explicit chiral symmetry breaking compared with 
the sum of the effects of explicit and dynamical 
chiral symmetry breaking. One anticipates that 
for light-quarks this ratio must vanish because the 
magnitude of their constituent-mass owes primar¬ 
ily to dynamical chiral symmetry breaking, while 
for heavy-quarks it should approach one. The fig¬ 
ure confirms these expectations. 

4. Summary and Conclusion 

We presented arguments which indicate that 
using solely spacelike information about a 
Schwinger function, it may only be possible to 
extract properties of bound states with masses 
< 1.2 GeV. With perfect spacelike information it 
is possible to accurately reconstruct the pole con¬ 
tributions from states which satisfy this bound. 
However, such reconstruction methods provide no 


information beyond that which is already avail¬ 
able in studies that employ physical (timelike) 
bound state momenta. We illustrated this via 
the Dyson-Schwinger equations. 

In the context of numerical simulations of 
lattice-regularised QCD, we speculate that it may 
be possible to reach higher mass states by em¬ 
ploying lattice data to constrain the infrared be¬ 
haviour of DSE integral equation kernels and sub¬ 
sequently using the DSEs to provide information 
on Schwinger functions at timelike momenta. 

The cr term is one useful gauge of the impact of 
explicit chiral symmetry breaking on a hadron’s 
mass. Such information is important for the in¬ 
terpretation of measurements that indicate a spa- 
cial and/or temporal variation in Nature’s funda¬ 
mental parameters. We calculated and reported 
(7-terms for a range of hadrons: the ground state 
pion; the pion’s first radial excitation; a light 
scalar meson; the p and w; the nucleon and A; 
and the u, s, c and b constituent-quarks. 

The analysis of the scalar channel is interest¬ 
ing. It is consistent with a picture of the lightest 
scalar as a bound state of a dressed-quark and 
-antiquark combined with a considerable two-pion 
component, which reduces the quark-core mass 
by < 15% but dramatically alter the current- 
quark-mass-dependence of the 0 ++ pole position. 

The constituent-quark er-terms provide a tool 
that is useful for assessing the flavour-dependence 
of dynamical chiral symmetry breaking. 
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